Interference in disordered systems: A particle in a complex random landscape 
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We consider a particle in one dimension submitted to amplitude and phase disorder. It can be 
mapped onto the complex Burgers equation, and provides a toy model for problems with interplay of 
interferences and disorder, such as the NSS model of hopping conductivity in disordered insulators 
and the Chalker-Coddington model for the (spin) quantum Hall effect. The model has three distinct 
phases: (I) a high-temperature or weak disorder phase, (II) a pinned phase for strong amplitude 
disorder, and (III) a diffusive phase for strong phase disorder, but weak amplitude disorder. We 
compute analytically the renormalized disorder correlator, equivalent to the Burgers velocity- velocity 
correlator at long times. In phase III, it assumes a universal form. For strong phase disorder, 
interference leads to a logarithmic singularity, related to zeroes of the partition sum, or poles of the 
complex Burgers velocity field. These results are valuable in the search for the adequate field theory 
for higher-dimensional systems. 



I. INTRODUCTION 

Much progress has been accompUshed in the under- 
standing of the thermodynamics of classical disordered 
systems [2 [2]. Typically, the disorder is modeled by a 
random potential. At low temperature, the low lying lo- 
cal minima of the resulting, rough, energy landscape be- 
come metastable states, and dominate the partition sum 
of the system. The correlations of the random potential 
determine the statistics of these metastable states, and 
hence the physics of the model. In many cases, for exam- 
ple in elastic disordered systems, the scaling of observ- 
ables can be described by a family of T = fixed points 
of the RG flow (most notably random-field and random- 
bond), which yield characteristic, universal roughness ex- 
ponents and effective disorder correlators [3HS]. 

However, the picture is much less clear when quantum 
interference is important. In real time dynamics one must 
study a sum over Feynman paths, whose weights are com- 
plex random numbers. The dominant contributions may 
then be difficult to discern. 

To give a specific example, hopping conductivity of 
electrons in disordered insulators in the strongly local- 
ized regime is described by the Nguyen-Spivak-Shklovskii 
(NSS) model |^. The probability amplitude J(a, h) is the 
sum over interfering directed paths F from a to 5 [Tllllj 
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The conductivity between sites a and h (e.g. on a T,'^ lat- 
tice) is then given by ^(a, h) ~ | J(a, &)p. Each lattice site 
i contributes a random sign r^j — ±1 (or, more generally 
a complex phase rjj = e^^'). 

Another example is the Chalker-Coddington model 
[T^ for the quantum Hall (and spin quantum Hall) effect. 
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where the transmission matrix T, and the conductance 
g{a, b) ~ tr T(a, h)^T{a, b), between two contacts a and b 
is given by [131 H] '■ 
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The random variables S'(ij) on every bond (ij) are U{N) 
matrices, with = 1 for the charge quantum Hall effect 
and A = 2 for the spin quantum Hall effect. Here F are 
paths subject to some rules imposed at the vertices. 

In both models, one would like to understand the dom- 
inating contributions to the sum over paths with random 
weights, given by J(a, b) or r(a, b), respectively; we shall 
denote it by Z in the following. In contrast to the ther- 
modynamics of classical models, where all contributions 
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FIG. 1: (Color online) Phase portrait of the model. The 
horizontal axis is the strength of V and the vertical axis the 
strength of 9. The effective disorder correlators for points A 
(deep in the diffusive phase), B (deep in the pinned ph ase) 
C (infinitesimally small disorder) are analyzed in sections III 
|IV| and [y] respectively. 
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are positive, contributions between paths with different 
phases can now cancel. One is also interested in the ex- 
pected phase transitions, i.e. critical values for the ampli- 
tude and phase disorder above which interference effects 
become important. 

This is a complicated problem. In this article, we 
therefore consider a toy model motivated by the mod- 
els above, for which many computations can be done ex- 
plicitly. More precisely, we analyze a "partition sum" Z 
defined in one dimension, of the form 



da;e"'^L^^"=^+'^^"=""'^'J"'^^"'^ (3) 



Here, V{x) is a random potential and 9{x) a random 
phase, both with translationally invariant correlations, 
and /3 = 1/T the inverse temperature. This is a toy 
model defined in one dimension and thus a drastic sim- 
plification compared to both the NSS model ([!]) and the 
Chalker-Coddington model ^ (which are usually con- 
sidered in two dimensions and above). However, a sim- 
ilar toy model for a random potential without random 
phases reproduces many important physical features of 
realistic, higher-dimensional models. For example, the 
appearance of shocks and depinning are already present 
in this framework |15H17j . For purely imaginary disor- 
der, a nice experimental realization of the sum in Q in 
cold atom physics is discussed in Section II B Complex 
sums similar to ^ are also encountered in magnetization 
relaxation in random magnetic fields, see e.g. jl8J. 

The interplay between random phases 0{x) and a ran- 
dom potential V{x) similar to ^ was already studied, 
among others, in [I9l - f22| . The basic distinction of three 
phases, high-temperature phase I, frozen phase II, and 
strong-interference phase HI, was established by Derrida 
in [53], and we follow his conventions. The aim of this 
paper is to pursue a complementary approach to [23], 
based on the study of renormalized disorder correlation 
functions. The latter is defined due to the presence of 
the parabolic well centered at x = w in ^. The resuh- 
ing spatial structure exhibits non-trivial features such as, 
in some cases, discontinuous jumps (shocks) as w is var- 
ied. Furthermore the renormalized disorder correlator 
is the central object of the field theoretic treatment of 
disordered systems based on the functional renormaliza- 
tion group [TH [21] , and thus the results of the toy model 
will give hints for a treatment of more realistic, higher- 
dimensional systems. 

This article is organized as follows: In section |Tl| we 
give the theoretical framework for our treatment. The 
model ([3]) is related to a complex Burgers equation, with 
time t — m^^, which has generated interest in the math- 
ematics community recently |25| . Equal-time correla- 
tion functions of the Burgers velocity field are the renor- 
malized disorder correlation functions A(w — w') of our 
model, ln(Z(u;))9u,' ln(Z(i(;')). The precise definition 
is given in Sec. jllDj They encode physical properties of 
the system like the appearance of shocks. Their m — > 



(i.e. < — ^ 00 in the Burgers picture) asymptotics forms 
the basis for the following analysis of the various phases. 

We first discuss the strong-interference regime, V{x) — 
and sufficiently strong 9{x), in section HI This is the 
regime most directly related to the NSS model and the 
Chalker-Coddington model described above. Naively, 
one may think in analogy to the case of classical disor- 
dered systems where 0{x) = 0, that points of stationary 
phase take on the role of the local minima of the energy 
landscape, and dominate the partition sum. We will show 
that this is incorrect. Instead, fiuctuations of Z{w) along 
the entire system are important. In our analytical treat- 
ment using the replica formalism, this manifests itself as 
a pairing of replicas. We will see that there is a finite 
density of zeroes of Z{w) (as already observed in [23]) 
which manifests itself in a logarithmic singularity of the 
effective disorder correlator A{w — w') for w close to w' . 

In section |IV[ we consider the infiuence of random 
phases in the frozen regime (large j3V{x)), where only 
a few local minima of the random potential contribute 
to Z{w). In the /? — > (X) limit one finds sharp jumps be- 
tween these minima as w is varied. In the Burgers veloc- 
ity profile, these manifest themselves as shocks, and their 
statistics are known to be encoded in a linear cusp of the 
effective disorder correlator [THldB]- We then discuss how 
the form of these shocks is modified by the introduction 
of random phases. It turns out that the linear cusp of the 
effective disorder correlator again acquires a logarithmic 
singularity. This time, however, it is related to shocks 
between two minima where the phase angle difference is 
TT, i.e. Z{w) passes through 0. This phenomenon is de- 
pendent on the spatial structure and on the possibility 
to vary and hence was not observed in [23!. 

In section |Vj for completeness we briefly discuss the 
high-temperature phase. Here, fiuctuations of Z are 
small compared to its expectation value, and Z never be- 
comes zero. As a consequence, the effective disorder cor- 
relators are regular everywhere, indicating that no shocks 
or poles occur in the Burgers velocity field. 

At any finite system size L ~ ^, there are blurred 
cross-overs between these phases as shown in figure [T] 
They become sharp transitions in the thermodynamic 
limit if the variance of 9 and V is rescaled with the sys- 
tem size L as ~ 6*2 ~ In L [53] . Since we are interested 
in the behaviour of the disorder correlators deep inside 
each individual phase, we do not follow this path but in- 
stead choose the simpler scaling ^ 1 or ~ L. By doing 
this for V and 9 individually, we can shrink all phases but 
one in the phase diagram to points respectively lines, and 
discuss each phase individually. 

In conclusion, one significant physical result of our 
work is that the introduction of random phases has quite 
different effects depending on the real potential V{x). If 
V{x) is sufficiently strong so that the system is in the 
frozen phase, even weak random-phase disorder imme- 
diately introduces zeroes of Z{w) or turns the shocks of 
the real Burgers velocity profile into poles of the complex 
Burgers velocity profile. In the high-temperature phase. 
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where V{x) is weak, this does not happen for weak ran- 
dom phase disorder, and the effective disorder correlators 
remain analytic. 



II. PRELIMINARIES 
A. Definition of the model 

To completely define the model (|3| , one needs to spec- 
ify the joint distribution of the random potential V{x) 
and the random phases 6{x). For the purpose of this 
paper, we assume that V and are independent, and 
that the distribution of 9 is symmetric around 0. This 
is mostly for technical reasons (since this choice makes 
many observables real) and is certainly true e.g. for cen- 
tered Gaussian distributions. 

Typically, one chooses V{x) to be Gaussian with mean 
zero, V{x) = 0, and variance depending on the type of 
correlations. Here denotes averages over realizations 
of the disorder. In the absence of imaginary disorder, a 
short-range correlated V{x), i.e. V{x)V{x') = unless 
X = x', gives the so-called Kida model [171 EH El]- If 
one chooses V{x) to be long-range correlated as a ran- 
dom walk, [V{x) — V{x')]'^ ~ |a; — x'\, one obtains the 
Sinai model IH]. For d{x) we also assume translationally 
invariant correlations. 

For some computations, it is easier to regularise by a 
finite system size L, 
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The system size L can be related to the mass of the har- 
monic well through L ^ ^ . In the case of pure random- 
phase disorder, i.e. V{x) = 0, Q can be seen as a parti- 
tion sum of a particle in the real random potential 9{x) 
at imaginary inverse temperature i/3. 



B. Proposed measurement of Z in cold atoms 

A direct measurement of the partition sum as given 
in ([3| for the strong-interference phase, i.e. with ran- 
dom phases 0[x) but without a random potential V{x)^ 
is at least in principle possible in a cold-atom experi- 
ment: Prepare the system in the ground state of a weak 
harmonic well (with frequency uj), so that sX t — Q the 
wavefunction is 



1 

M^) :=^(x,t = 0) - (^^y 



(5) 



Then, switch off the harmonic well, and instead switch on 
a random potential 0(x). In situations where the kinetic 
term in the Hamiltonian is negligible, such as <C 1, 
or a large mass M , the time evolution is approximately 



given by e i.e 

ip{x,t) 



/Muj 



e 2R 



-j-e{x)t 



(6) 



Switching the potential back to the harmonic well and 
measuring the overlap with the ground state gives 



dx ( 



(7) 



This is exactly of the form of Z{w) in (|3|. Although 
the overlap ('0o|V'(^)) cannot be measured directly, the 
occupation probability of the ground state, given by 
I (■(/'o I "0(0)1^1 could in principle be measured, providing a 
direct measurement of |Z(0)p in the strong-interference 
regime. An example of a related experiment is given in 



C. Complex Burgers equation 

Another application of the toy model (|3| is to the com- 
plex Burgers equation. With the mapping t :— m~^, 
one obtains from ([s]) the "equation of motion" or "renor- 
malisation group equation" {m^ being interpreted as the 
infrared cutoff) for Z, 



dtZ{w,t)^^dlZiw,t) . 



(8) 



We have added the argument t for clarity. The initial 
condition at t — equivalent to = cx) is 



Z(w,0) = e-'3^(«')-*«("') . 

Using the Cole-Hopf transformation h{w,t) 
~T In Z{w,t), we obtain the KPZ-equation 



dth{w,t) = ^dlh{w,t) 



(9) 



(10) 



Taking one spatial derivative one arrives at the Burgers 
equation for the velocity u{w,t) := du]h(w,t): 



T 

dtu{w,t) = -dlu{w,t) 



ud,ji]U{w, t) 



(11) 



Without random phases {6{x) — 0), u{w,t) is real and 
(11) is the well-studied real Burgers equation. It has 



been used, among others, to describe the formation of 
large-scale structures in cosmology (the so-called adhe- 
sion model) and in compressible fluid dynamics (for a 
review see [30J). When random phases are included, u be- 
comes complex. The resulting complex Burgers equation 
has very surprising applications e.g. to Lozenge tilings of 
polygons [5S]. It has also been used to obtain further in- 
formation on the real Burgers equation through analytic 
continuation to the complex plane and the so-called pole 
expansion |31H34I . In the following, we study the small- 
m properties of ([3]) , or equivalently the large-t properties 
of the complex Burgers equation ( 11 ). 
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D. Effective disorder correlators 

The main observables on which we base our analysis 
are the so-cahed effective disorder correlators, which we 
define now. For each reahzation of the random potential 
V{x) and the random phases 0{x), we first define the 
"free energy" or the "effective potential" by 



pV{w) + i9{w) := I3h{w) = - In Z{w) 



(12) 



Note that V{w) is always unique, but 9{w) is only de- 
fined modulo 27r. We will thus focus on 9'{'w) which is 
unambiguous. This is also the reason why it is prefer- 



able to consider the Burgers equation (11 1 instead of the 
equation ( 10 ) for the potential. 



The effective disorder correlators for the potential and 
the phase are then defined by 



Av{wi-W2) ■.^V'{wi)V'{w2) 



Ag{wi-W2) ■.= 9'{wi)e'{w2) . 



(13) 



The cross-correlator V'{wi)9'{w2) vanishes since V{x) 
and 9{x) are independent and due to the symmetry 
9^—9. In more general situations, this may not hold. 

The correlators defined above have a nice represen- 
tation as correlation functions. Define the normalized 
expectation value of an observable O for a given w as 



{0{x))^ :^ 



1 //3m2 
Z{w) V 



(14) 

-'^'■'^^0{x) . 



By definition (1)^ = 1. Taking a derivative of (12) yields 
pV'{w) + i9'{w) = /3m2 {w - x)^ 



(15) 



This gives two simple relations for Ay and A^: 

4 



Azz(wi - W2) := m {x - wi)^^ {x - wa)^, 

= Ay(w;i - W2) ~ (3~'^Ao{wi - W2) 
Azz' {wx - W2) := m^i^x - w^)^^ {x - W2)l^ 
= Av{wi-w2)+l3-^Ag{w 



(16) 



(17) 

■W2) 



In terms of the complex Burgers equation ( 11 ), the effec- 



tive disorder correlators have the intuitive interpretation 
of equal-time velocity correlation functions: 



Azziwi - W2) = u{wi,t)u{w2,t) 
Azz'iwi - W2) = u{wi,t)u*{w2,t) 

We now proceed to compute Azz and Azz* in each of 
the three phases, and discuss their implications on the 
physics. 



III. THE STRONG-INTERFERENCE PHASE 
(PHASE III) 

The strong-interference phase has first been discussed 
in the context of directed paths with random complex 
weights in [13 HH Ei] and later for the random-energy 
model at complex temperature |23j . In this phase, the 
average of Z is essentially zero (or at least subdominant) 
due to strong interference, and Z is dominated by fluc- 
tuations. The whole system contributes to the partition 
sum, in contrast to the case of a real random potential, 
where it is dominated by a few points with exceptionally 
large moduli. 

In a replica formalism, this is reflected by a pairing of 
the replicas, as already observed for the NSS model in 
|35j . For the two-dimensional model discussed there, an 
entropic attraction between replica pairs arises at cross- 
ings of four or more replicas due to the spatial structure. 
In our one-dimensional model, the resulting replica pairs 
will turn out to be essentially non-interacting and spread 
out over the whole system. 

We will analyze this phase by setting V{x) — in 
([3]) and consider the small-m limit. We shall show that: 
(TJ this phase is characterized by Z{w) being a Gaus- 
sian stochastic process in the complex plane with w as 
the time variable; (ii) its two-time correlation function is 
universal and given by 



Z{w)Z{w') - e — r 



8m^ ( Isl 



(18) 



From this, the effective disorder correlators defined above 
can be computed. We shall see that Ay and Aq exhibit a 
logarithmic singularity around zero, describing the statis- 
tics of zeroes of Z. In contrast, Azz — Ay — /J^'^Ag 
remains regular around zero. 

We then consider two explicit examples where the 
random phase disorder is sufficiently strong to observe 
this phase. Example 1 will be a model with Brown- 
ian imaginary disorder, i.e. long-range correlated phases 
\9{x) — 9(x'y^ ^ \x — x'\. Example 2 will be a model 
with short-range correlated phases uniformly distributed 



A. Characterization of phase HI and probability 
distribution of Z 

We set V{x) = in (|3| and consider imaginary dis- 
order. There is a large class of processes 9(x) such that 
in the limit m — >■ the distribution of Z{w) tends to a 
complex Gaussian variable due to a central limit theorem 
(CLT). To understand qualitatively why, let us think of 
Z(w) as a discrete sum Z{w) ~ J2j=i where each 
Zj — e*^J is a random variable inside the unit disk and 
L ^ 1/m. The usual statement of the CLT shows that 
uncorrelated variables Zj belong to this class, (this is ap- 
plied e.g. in example 2, section HIE I. In the more general 
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case of correlated Zj a CLT also holds under the assump- 
tion that the correlations of the Zj decay fast enough. 
A precise mathematical statement of the necessary and 
sufhcient conditions is possible using a so-called strong 
mixing condition (sec 36: ,38j). 

More qualitatively, we require the conditions that 



tion group" equation (Is]) describes how Z{w) evolves un- 
der changes of t = m . This implies a similar equation 
for the 2-point function f{w — w') :~ Z{w)Z*{w'), 



dtft{w)^TdlMw) 



(20) 



q± = 



(19) 



are finite and that a similar condition for the integral of 
the fourth cumulant holds. Since we assumed that 0{x) is 
symmetrically distributed around 0, the q± are real. Note 
that the fact that Z{w) is bounded for any realization of 
0{x) and any w by |Z(z/;)| < 1 distinguishes this problem 
from the real potential case, where the CLT does not hold 
in general. 

Thus, from now on we consider the case where the CLT 
holds and in the limit m -> the distribution of Z(w) 
tends to a complex Gaussian variable. This happens in 
what we call phases I and IIL In these phases, the distri- 
bution of Z{w) is thus determined by its mean Z(w) and 
its covariance matrix, consisting of 3 entries Z(w)Z{w), 
Z{w)Z* (w), and Z*(w)Z*{w). A similar reasoning ap- 
plies to the joint distribution of ZivS) and Z{w''). 

The key difference between the strong-interference 
phase III and the high-temperature phase I is the scal- 
ing of these moments: If the mean of Z as a function 

2 

of m decreases faster than the fluctuations, Z(w) <C 

(Z(w) — Z{w))^ ^ m as TO — >■ 0, we obtain the strong- 
interference, fluctuation-dominated phase III. If, on the 
other hand, the mean decreases slower than the fluc- 

0, we obtain the high- 



tuations, Z{w) 3> to as to 

temperature phase I. 

In the example in s ection |IIID| w e will take 9 to be 

long-range correlated, [9{x) — 0(y)]^ ^ \x ~ y\, whence 

2 _ , 

we will see that Z{w) ~ e ^ to. On the other 

hand, in section [ill E[ we will consider an example where 

rotational symmetry enforces Z{w) = 0. In both cases, 

we verify the general results and assumptions presented 

here. 



B. The second moments 

The second moments of the complex process Z{w) take 
a general form which we derive now. The "renormaliza- 



Here, we added the index t to make the dependence of 
Z(w) on the parameter to and hence the dependence of 
f{w) on the parameter t explicit. The general solution of 



(201 in terms of the initial condition fo{w) = e'(^(o)-^('")) 



IS 



AnTt 



e f^(^wo) dwo . (21) 



Since we assumed q- to be finite, see Eq. (19), the solu- 



tion (21 ) tends to a Gaussian scaling form as t — > oo 



ftiw) 



AnTt 



(22) 



Note that this assumption is violated in phase I, where 
the mean Z{w) contributes a constant to ft{w) even 
for large t. It is also violated in phase II, where 
I^oo fo{wo)dwa diverges. 

Going back to the original variables to and /3, the 
asymptotic form of (21 1 in phase III is 



fm{w 



^"^/(^ = my^{w - w')) (23) 
. (24) 



Exactly the same re asoning goes through for the second 
moment Z{w)Z{w') with q_ replaced by q+. 

The scaling in (l23|) reflects the fluctuation-driven char- 



acter of phase III: If the mean Z{'w) were not subdom- 
inant, for large system sizes L ~ m~^, as compared to 
the fluctuations, / in (23) would be 0{LP) instead of 



0{L ^), and not tend to zero for large argument. 

To summarize, in the strong-interference phase III, as 
TO — 0, the partition function Z{w) tends to a Gaussian 
process with mean and correlation matrix: 



/ Z{w)Z{w) Z{w)Z*{w) Z{w)Z(w') Z{w)Z*(w') \ 
Z*{w)Z{w) Z*{w)Z*{w) Z{w)Z{w') Z{w)Z*{w') 
Z{w')Z{w) Z{w')Z*{w) Z{w')Z{w') Z{w')Z*{w') 
\Z*{w')Z{w) Z*{w')Z*{w) Z*{w')Z{w') Z*{w')Z*{w')j 



( q+ q- q+f{w) q-f{w) 
q- q+ q-f{w) q+.f{w) 

■ q+f{w) q-f{w) q+ q- 
\q-fiw) q+fiw) (7_ q+ 



(25) 
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With this, we have completely characterized the m — 
asymptotics of Z^w) in the strong-interference phase 
III as a Gaussian stochastic process with the second mo- 



ment given by (23 1. In sections HID and HIE we shall 



explicitely check the asymptotic form in ( 24 ) and obtain 
the non-universal constant q± in (24 1. 



The disorder correlators 



Plugging in the scaling form /{w) = e * 
(29), we obtain the disorder correlators 



into (28) and 



Av{w) = 
Equivalently, 



m 



e^"^^ - 1 



-^I^ln(l 



(30) 



Having discussed the probability distribution of Z{w), 
we now turn to computing the effective disorder corre- 
lators. For simplicity, we restrict ourselves to the case 
when the limiting Gaussian distribution for Z(w) is ro- 
tationally symmetric, i.e. only depends on the modulus 
|Z(u))|. In the covariance matrix (25 ), this means g+ = 0. 



The joint probability distribution for two partition sums 
Z{wi) = ai+ihi^ and Z{w2) = 02 + i 62 is then given by 



^'(ai, 61, 02, 62) = 



with x = ( oi 61 02 62 



1 

47rVdetS 



(26) 



1 f{w) 
1 f{w) 
2 I f{w) 1 
f{w) 1 



(27) 



and c — mq^ J ^ . We will see later that the disorder 
correlator does not depend on c. 

To compute the effective disorder correlators, let us re- 



consider their definition ( 13 ). Since 6 is the angular vari- 



able of a two-dimensional Gaussian stochastic process, we 
can apply the results of [17 . There, the two-time correla- 
tion function for the angular "velocity" 9{'w) := dw0(w) 
of planar Brownian motion is (cf. [39], formula 17): 

§{w)kw') = -^[d^d^>lnf{w)]ln{l - f{wf) . (28) 



The two-point correlator of the phase (instead of its ve- 
locity) can then be written as a double integral of (28), 
but no explicit expression is known. 

The two-point correlator of ln|Z| = (3V is obtained 



from the explicit form ( 26 ) for the two-time probability 
distribution as 

P'^V{w)V{w') = \n \ Z{w)\ In \Z{w')\ 



-00 47rVdctB 



e 2 



In lai -I- i &i| In |a2 -I- i 62I 



with B given by (27). This integral can be computed 



exactly (7E denotes Euler's constant): 



l3^V{w)V{w') 



1 r 

4 



(7E-lnc)VLi2(/(u;)2) (29) 



Azz{w) = - 



4/3 65-' - 1 



Azz'{w) 



m 



w 



g2 IB 



-+21n( 



(31) 
(32) 



Observe that Azz is smooth around 0, whereas Azz* 
has a logarithmic singularity at w = 0. Note that all 
correlators are expressed in terms of the rescaled variable 



w defined in ( 23 ) . 



The above expressions for the correlators (which are 
also the two-point equal-time velocity correlators for the 
complex Burgers equation) are universal and generally 
valid in phase III, under the assumption of rotationally 
invariant disorder. The more general case can be handled 
by similar methods but is not studied here. These results 
were obtained using the CLT assumption. 

We now study two specific models where we can com- 
pute the general moments (beyond the second one) using 
the replica method, and check that they are consistent 
with the above reasoning. As an additional check we 
also compute numerically the correlators. 



D. Example 1: Imaginary Brownian disorder 

Consider pure random-phase disorder, V{x) — 0, and 
take 0{x) to be a continuous random walk, i.e. a Gaussian 
stochastic process satisfying 



(33) 



Thus in the finite length rcgularization, the partition sum 
^ is the Sinai model at an imaginary temperature. We 
measure numerically the effective disorder correlators, us- 
ing relations (16) and (17). The results are compared in 
figure [2] against the analytic computation in the previous 
section. We observe good agreement. 



1. Second moment 

Here we show explicitely the validity of the scaling ar- 
gument given in section [ill B| for this model. 
Using formula ([3|, the second moment is given by 



Z{w)Z*{w') 



(34) 



2tt 



dx dy e 



-^a\x-y\-l3!!^[(x-wf + {y-w')'' 



A(w)//?m- 
2.0 




FIG. 2: Effective for ce-force correlators for the long-range 
model defined by (33 1. Red line: Azz from 131 1, blue line: 



Azz* from (321. Dots: correlators obtained from numerical 
simulations using ( |16[ ) and | |17[ ) for a = 1,13 = 10 and m = 
0.05 (blue and purple), m = 0.1 (yellow and green). 



This integral ( 34 ) can be computed exactly by using the 



center-of-mass variable s and the "pair separation" vari- 
able t, defined as 

x + y 

t :— X — y . (35) 
In these variables, the s and ^-integrals decouple, 
/3m2 ^ 



Z{w)Z*{w') = 



2tt 



(itQ-\''\t\+l3'^{w-w')t-l3^t^ 



(36) 

Both integrals can be performed analytically. In terms 
of the rescaled variables 



w := my/p{w — w') 
a 



(37) 



my/13 



the second moment ( 36 1 is given by 



Z{w)Z*{w') 



r cr-w 

e 1- Erfc 

2 

I a + w 

+e 4 Erfc 



(38) 



Taking the m — >■ limit at fixed (3 (i.e. the limit a — > oo), 
the second moment ( 38 1 approaches the scaling form 



Z{w)Z*{w') 



(39) 



This is exactly the scaling form obtained in (24), and 



2. Higher moments 



Let us now look at higher moments of Z and Z* , given 



by: 



Z[wi)...Z{w,,)Z*{w[)...Z*{w'^) 

oo oo 

= / dxi • • • Axn I dj/i • • • dyn 



xe^'^T^ Er=i[(2:i-"'x)^+(2/.-«'i)^] (40) 

An exact calculation does not seem feasible, but the 
asymptotic behaviour is understood as follows: In the 
limit (7 — > oo, the exponent in (40 1 will have a sharp max- 



imum at configurations where the Xi and yj are paired, 
i.e. close to each others. We now consider configurations 
which are close to such a pairing, where w.l.o.g. Xi is 
paired to y7r(i) with some permutation tt. Similar to (35) 



we introduce center-of-mass and separation coordinates 
Si and ti for each pair and rewrite the mass terms as in 



(36) 



The ^i-integrals have complicated boundaries, which 
yield terms decaying as e~°"^ with various functions a > 
0. Hence, these terms can be neglected in the limit m 
0, and the s and t-integrals decouple again: 



Z{wi)...Z{wn)Z*{w'{)...Z*{w'^) 

n 

= ^ Z{wi)Z*{w'^^-^) + higher orders in m (41) 
In particular, we get 



[Z{w)Z*{w')f =nl Z{w)Z*{w') 

+ higher orders in m 



(42) 



A more rigorous justification that this is the leading term 
in an expansion in orders of m is given in appendix [A| by 
considering the moments of the partition sum in a finite 
system (|4]). 

Correspondingly, the leading term for the moments 
[Ziw)]^ \Z* {w')]™" for m 7^ n is zero in the strong- 
disorder limit, since then the replicas cannot be paired. 
Stated differently, the phase of Z is random, and hence 
only moments invariant under the rotation Z — )> e^'^Z 
are nonzero. Dropping the higher-order terms in (41), 



gives a non-trivial check for the validity of that argument. 



we obtain exactly the moments of a complex Gaussian 
variable. This supports the general claim made in sec- 
tion |III A[ and shows that this model is indeed in the 
strong-interference phase III. 

The fact that configurations with unpaired replicas are 
subdominant shows that fluctuations of Z dominate over 
the average. Intuitively, this happens since for ct 3> 1 
the phase of the integrand in the expression ([s]) grows 
beyond 2tt on a scale much smaller than the width ^ 
of the parabolic well. Hence, (|3| is essentially a sum of 
many random complex numbers with mean 0. 
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For n = 1, evaluating the second moment gives: 



3 4 5 



FIG. 3: Effective forc e-force correlators for t he short-range 
model defined by (43 1. Red line: Azz from 131 1, blue line: 
Azz* from (321. Dots: correlators obtained from numerical 
simulations using ( |16[ ) and | |17[ ) for a = 1,13 = 10 and m = 
0.05 (blue and purple), m = 0.1 (yellow and green). 



This is the same behaviour as in "phase III" discussed 
by Cook and Derrida 19J and by Derrida [23]. Since 
our potential is long-range correlated, 6{x)'^ ~ x instead 
of the short-range correlated potential 9{xY ^ 1 used 
in [in] , the complex phase of the integrand in ^ grows 
much faster in our model. Hence, we do not observe 
"phase I" for high temperatures (/3 < /3c) as in [T^ I23j . 
but only the fluctuation-dominated "phase III" . 

In the following, we shall show that similar results hold 
in a model with uniformly distributed 9(x). 



E. Example 2: A short-range correlated model 
with uniformly distributed angles 



Z{w)Z*{w') 



27r 



^--I3^[{x~wf + {x-w'f] 



(44) 



Note that this is again in agreement with the general 



scaling argument given in section III B[ 

For the higher moments [Z{w)Z* (w')]", the only terms 
contributing are those where the 2n replica form n pairs. 
When at least two pairs are at the same position, it is not 
clear which replicas are pairs, leading to double-counting. 
However these contributions are subdominant and vanish 
with a relative factor of at least m?. Thus, the dominant 
term for m — > is 



[z{w)z*{w'y\ 



= T7 ' 



[3m? 



-li^[{x~w)^ + (x-^w'f]^^ 



Z(w)Z* (w') + higher orders in m 



Analogously to the derivation of (42) this formulas can 



be generalized to moments of Z with different positions. 

Again, we observe the same behaviour of the moments 
as for a complex Gaussian. In total, in the limit of large 
m, we recover the same phase HI results as in the long- 



range correlated model in section HID and confirm the 
validity of the general arguments in the beginning of this 
section once more. 



IV. THE FROZEN PHASE (PHASE II) 



Our second example is a model, where the potential 
9{x) in ^ is short-range correlated and uniformly dis- 
tributed. To be more precise. 



p{e{x)) = 



27r 



(43) 



while 9{x) and 0{x') are uncorrelated for x ^ x' . 

As can be seen on figure |3| a numerical simulation 
yields disorder correlators which compare well to the gen- 
eral results obtained above. As for the first example, we 
shall compute moments of Z to elucidate the physics. 

Invariance of the distribution of Q(x) under a phase 
shift, Q(x) — ?> B(x) + 0, implies invariance of the distri- 
bution of Z under Z — >■ Z e^"^ . Hence, the only nonzero 
moments are of the form |Z(w)p"'. 



For large /3 and sufficiently strong potential V{x), the 
modulus of the integrand in ([3]) has a very broad distri- 
bution. It is well-known that the partition sum (in the 
absence of the harmonic well) is then dominated by a 
few points, the minima of Vix). This so-called frozen 
phase has been extensively studied in the absence of ran- 
dom phases by a variety of methods (replica symmetry 
breaking [27J, functional renormalization group [17] and 
rigorous mathematical analysis [40J) 

Distributions of V where a frozen phase occurs in the 
model ([3]) in absence of complex phases include: 



• Long-range correlated random potentials V{x), i.e. 

1. This is known as the Sinai 



V{x)V{x') =a\x-x\ 
model, which describes the diffusion of a random 
walker in a ID random static force field [551 EI] • 

• Short-range correlated random potentials V{x), i.e. 
V{x)V{x') = ad{x — x'), where the amplitude is 
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rescaled logarithmically with the system size, or m: 
a — Inm. Freezing occurs below some critical 
temperature, j3 > /3c, analogously to the random 
energy model [ID]. 

Among the most interesting features of the frozen phase 
is the appearance of jumps between distant minima of 
V{x) as the position w of the harmonic well in ^ is 
varied [121 fi^Hl^ . These correspond to shocks 30 under 
the mapping to the Burgers equation discussed in section 
|II C[ In the following, we will discuss how their structure 
is changed upon introduction of random complex phases 
0{x), following the standard treatment [T71 |3S] for the 
case without random phases. 



A. Complex shocks 

Let us first consider a fixed realization of the random 
potential V{x) and the random phases 9{x). For almost 
all w, the real part of the exponent, V{x) + ^{x — w)'^ 
has, as a function of a;, a single minimum at some value 
X = Xm{'w). Then, in the low-temperature limit (i.e. 
/3 ^ oo) 



Z{w) 



and hence 



V{w) = V{x,n) + ^(a;m - w f 
9{w) = e{x^) 



(46) 
(47) 



The function Xm{w) is constant over some range of w, but 
then jumps to a different value at w — w* . Denoting the 
two solutions at w* by xi and X2, the necessary condition 
for a jump is 



^(^i) + ^(a^i - w*? = V{X2) + —{X2 - w*f (48) 

In terms of the effective potential V , two parabolic sec- 
tions given by (46) (with w = wi and 'W2^ respectively) 
meet at w* with a linear cusp. The first derivative, 
V'{w), has a discontinuity at w*. 



So far, this is the same picture as has been established 
for purely real disorder long ago in the context of the 
Burgers equation [T3 EHl ISj- There, the appearance 
of the shocks is succintly encoded [TH [II] in the effective 
force- force correlator A (it;), which extends to the broader 
context of interfaces in random media. It has been com- 
puted and tested both numerically [17| EH] and exper- 
imentally |49j. It encodes the statistics of the shocks 
through a linear cusp aX w — 0. At finite temperature /3, 
the shock is smoothened in the so-called thermal bound- 
ary layer which extends on a scale w ~ T = /3~^ [50ll5T] . 

The additional random phase 9{x) will in general be 
different at Xi and X2. We now show that this is reflected 
in the profile of V{w) and 9{w) for w close to a shock. 



This modifies the form of the disorder correlator A(it;) 
near w = 0, more specifically in the thermal boundary 
layer region w ^ T where we will obtain its precise form. 
We find that it adds a logarithmic singularity which de- 
pends on the statistics of the phase jumps. 



The shock profile 



Let us assume a two-well picture, i.e. approximate 

Z{w) by 

Z{w) = p-l3[Vi + ^{xi-w)^]-ie, ^ ^-i3[V2 + ^{x2-wf]-te2 

(49) 



The effective potential ( 12 ) can be written in terms of 



the jump size s := I3m?{x2 — xi), the phase difference, 
(/) :— 92 — 9i, and w* , solution of (tlsl): 



9\w) 
-V'iw) 



sin(0) 



2 cos(0) -|- cosh {s[w — w*]) 
s sinh (s[w — w*]) 
2/3 cos(0) -t- cosh {s[w — w*] 

+ ^{^1 +X2- 2w) 



(50) 
(51) 



(45) Some examples of shock profiles for various values of the 



parameters are shown in figure IV A 1 Note that as — >■ 
±7r, a pole appears in V'{w) at w — w* which is the real 
part of the Burgers velocity. 

To obtain the disorder correlator Ag, we need to aver- 
age 9' {wi)9' {W2) over the disorder. Assume a small uni- 
form density po of shocks, and average over w* with the 
1 

measure po /^'!L dw*. Since 9'{w) decays rapidly as w* 

is increased, we can safely extend the integration limits 
to ±cxD, allowing to compute the integral over the shock 
position w* analytically: 

/7r />oc 
A4> dsP{cb,s)f{<j>,s) (52) 



/((/>, s) = ssm^{(t>)- 



)COt( 



^ coth 



COS 2(p — cosh sw 



(53) 



where w = wi — W2- We denote by P(4>,s) the joint 
distribution of the jump sizes s and the phase jumps (p. 
Remarkably, Ay{w) can also be calculated, by consider- 
ing the difference Ay{w) — Ay,p^Q{w), where Ay,p^q(w) 
is the correlator of the problem without the imaginary 
disorder, 9{x) — 0: 



Av{w) = Av,4,=oiw) + r^Ag{w) . 



(54) 



Thus the correlator Azz = Ay — (3~'^Ag is unchanged 
by the complex phases. Observe that the integrand in 
(52) becomes singular for w = and (f) — ±7r. In the 



following examples (sections IV B and IV C ) we shall see 
that this singularity yields a logarithmic singularity in 
Ay£ around zero. Its coefficient will be seen in section 
IV C to be proportional to P{(j) = ±7r). 
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0'(w) 
3.oi[ 

2 
2: 
1.: 
1 




(a) 



FIG. 4: Shock profiles from (51 1 for P — 5, w* = 0, xi + xi 
(dot-dashed), — ^-k (sohd). 




= f (dashed), </!) = fyr 



We have now discussed the effective disorder correla- 
tors Ag and Ay in a two-weU approximation in a general 
situation. So far, we did not make specific assumptions 
on the distribution and the correlations of the disorder. 



These enter the final result (52 1 through the joint distri- 
bution of the jump sizes s and the phase jumps 4>- Now, 
we will specialize to examples of particular interest. 



B. Example 1: Uniformly distributed random 
phases in a short-range potential 

Our first example is B{x) uniformly distributed in 
[— TT, tt] and uncorrelated from the spatial dependence x^, 
i.e. P(0, s) — ^P{s)- This allows to perform the </> inte- 
gral in ( 52 ) analytically: 



Ae(w) =Po / dsP(s) 



- 1 



ln(l 



(55) 

To take the limit /3 — >■ oo, we write s — (im?{x2 — xi) — 
/3mfis, where /i is the jump-size scale and the distribution 
P{s) is known as the Kida distribution [171 HZ] , 



P{s) 



-se * 



(56) 



Observe that the scaling is different compared to the cor- 
relator in the strong-interference phase III: The argu- 
ment of the scaling function is now w = Pmiiw instead 
oi w ~ in (30). 

For small x, (58) has the asymptotic form 



^e(a;) = 7(7: 



21n: 



Oix) 



(59) 



This logarithmic singularity arises from the — ±7r limit 
of the integral ( 52 ) , and is hence caused by shocks where 
Z(w* 



0. 



The integral ( 58 ) can be computed numerically and 
compared to simulations. We obtain a very good agree- 
ment with our numerical results (see figure [s]) for various 
values of m and /3, providing a non-trivial check for the 
scaling in (58). The scale /i is fitted as /i = 0.58, inde- 



pendent of m or /3 in the considered range 



Example 2: Wrapped gaussian distribution in a 
short-range potential 



The scale is related to the density of shocks po through 

m 

i = Po {x2 - xi) ^ po— dsP(s)s = po — Vtt . (57) 



We thus obtain the scaling form 
Ae(w) = (3ni^Ag {w — (im^w) 



(58) 



ds 



ln(l- 

1 ^ 



It is interesting to consider an example where the dis- 
tribution for (j) is non-uniform. We take again the phase 
angle 9{x) to be uncorrelated at different points. At each 
point, we assume the distribution of G [— 7r;7r] to be a 



^ ActuaUy, for the short-range random potential on a discrete lat- 
tice considered here, contains corrections, which are scaling 
logarithmically with m; see 1171 for more details. If we were to 
perform the simulations with m varying over several decades, ^ 
would have to be adjusted correspondingly. 
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FIG. 5; (Color online) Correlator in phase II, with 9{x) uni- 
form in [— TT, tt]. Dots: simulation results (blue: m — 0.01, 
P = 10, red: m = 0.01, /3 = 20, green: m = 0.05, /3 = 20, 
black: m = 0.05, /3 = 40, orange: m = 0.05, (3 — 60), rescaled 



as in (58 1 with n = 0.58 (for all curves). Blue line: (58 1, red 
line: asymptotics (59 1. 



2.0 a 



1.5 



1.0 



0.5 



0.0 




FIG. 



in (601 



6: Correlator in phase 2, with 9{x) wrapped gaussian as 
Dots: simulations (From top to bottom the variance 
-- 2, a ^ 1.6, a = 1.2, a ^ I, a = 0.8, 
0.01, and /3 — 120. Correspondingly, the 



decreases as: a ■ 
the mass is m = 
probability of a jump through zero decreases as P{(j} — ±7r) = 
0.15,0.13,0.08,0.05,0.01.), rescaled as in ^ with ^ = 0.59 
(for all curves 



Lines: numerical integration of (52 1 



wrapped gaussian distribution with variance 



1 1 o2 



(60) 



2ni 



■d denotes the Jacobi theta function. Note that P{9) is 
periodic, P{9 + 2n) = P{6). From ([60]), the distribution 
of phase jumps <j) — 62 — di is 



P{(t>) = I P{9)P{0 + (j3)d0 



(61) 



For the random potential, we still assume a short-range 
random potential as in section [IVB| Th e distribution of 
jump sizes s is thus still given by ( [56^ . This allows to 
obtain the full disorder correlator Ag by computing the 
integral (52) numerically. The results in figure [6] compare 
well to numerical simulations. 

One again observes a distinctive logarithmic singular- 
ity at w = 0. This arises from the (p = ±7r limit of the 
integral (52). More precisely. 



d(/)F(</>) sin^(0) 
-7rP(, 



(j) cot (j> - 



f coth 



COS 20 — cosh sw 
±7r) \nw + 0{w°) 



(62) 



The integral over s is normalized since 



r 

JO 



As = 1, and hence 



Kg{x) = -■nP{(j) = ±7r) Ina; + 0{x^) . 



(63) 



For a uniform distribution of 9, P{(j) = ±7r) = ^ and 
we recover the In part of the result (59). The constant 



In general, the coefficient of the logarithmic singularity 
at = is proportional to the probability of phase jumps 
by an angle of = ±7r. Thus intuitively this singularity is 
caused by shocks between minima of V{x), where Z(w* — 
0+) changes to Z{w* + 0+) = -kZ{w* - 0+) where k is 
a positive number. Note that at any temperature T > 0, 
i.e. /? < 00 the function Z{w) is smooth, thus passes 
through zero in our two-well approximation. This means 
that the Burgers velocity profile has a pole. Observe that 
according to (63), the logarithmic singularity is present 



as soon as there is a finite probability of jumps with an 
angle oi cj) — ±tt, however small it may be. This shows 
that in our model, there is no "sign phase transition" in 
the frozen, or pinned phase. This is in agreement with a 
recent result for a higher-dimensional model [52] where 
only 9 = and 9 = tt, i.e. plus and minus signs were 
considered. 

It is straightforward to repeat the analysis above with 
long-range correlated random potentials V{x). For ex- 
ample, in the case of the Sinai model explicit expressions 
for the probability distribution (56) for the jump sizes s 



coefficient of order w is harder to obtain. 



are known (see [ITj) but lead to complicated integrals. 

Note that the two-well model is only expected to be 
valid asymptotically for /3 — >■ cxi. At low but non- 
zero temperature, we expect subdominant contributions 
from higher-lying minima which may provide additional 
rounding of the singularities discussed above. 

In the following, we shall see that the behaviour in the 
high-temperature phase is quite different. 
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V. THE HIGH-TEMPERATURE PHASE 
(PHASE I) 



A(w) 

0.15 r 



For completeness, we also discuss the disorder correla- 
tors in the high-temperature phase. In this phase Z is 
dominated by the average and fluctuations are sub- 
dominant. As a consequence, for example the quenched 
average of the free energy is equal to the annealed average 
of the free energy. 

In our one-dimensional model, this phase occurs for 
sufficiently weak random potentials (for example, short- 
range correlated Vix) below a critical value of /?, which 
increases with system size) and sufficiently weak random- 
phase disorder (for example, short-range correlated Q(x) 
with finite variance, e.g. a wrapped Gaussian distribu- 
tion) . 

To compute the leading-order term for the correlators, 
let us take the example of short-range real and imaginary 
disorder, with 



V{x) = 0{x) = 
V{x)V{x') = av6{x - x') 



{x)9{x') = (Jg5{x - x') 



(64) 
(65) 
(66) 



For small ay and ag, we can expand the partition sum 
in powers of V and 0: 



Z[w) 



2 /-oo 



(67) 

dx [1 - (3V{x) - i9{x) + ...] e-^T^(^''")' 



The leading order for the effective potential thus be- 
comes: 



Viw) = \l^ r y{x)e-<^'^^''-^^^ Ax (68) 



27r 



CO 

2 /-oo 



6'(a;)e-'3'V(^-«')'da; (69) 



From this, we obtain the leading order for the disorder 
correlators in the high-temperature phase: 



Ay(w) — ay 



(2 - w^) e"T" 



(70) 



. , ^ (Bm 1^ , 9N ™- , , 

^e(w)^ae ^' ' {2^w^)e-^ (71) 
oy TT 

Here, w — m^/]3w. 

Another way to understand these correlators is through 
the so-called exact renormalization group equations fol- 
lowing [13. From ([12]) and (|3|, we obtain a flow equation 
of the form: 



mdmViw) = -^dlViw) - ^ ( d^V{w) ) (72) 



For the correlator R{w — w') := V{w)V{w'), this gives: 
- mdmRiw) = -^diR{w) + ^^iio(0, 0, w) (73) 



0.10 



0.05 



0.00 



-0.05 



FIG. 7; (Color online) Renormalized disorder correlators in 
the high-temperature phase. Dots: Simulation results for /? = 
0.1, av = 1, ae = 0.25, m = 0.5. Blue (red) line: Ay (As), 



as obtained from (701. No fit parameter 



Here, S{'Wi,W2,W3) :— V{wi)V{w2)V{w3) is the third 
cumulant and the subscript S'no indicates derivatives 
with respect to the first two arguments (notations as in 
jl7j). Without the non-linear term, equation (73) is the 



same as ( 20 ) , solved by ( 22 ) with initial conditions ( 65 1 , 



I.e. 



R[w)=ay\ — — e ^ 
47r 



(74) 



The feeding term for 5* is of order RR ~ /3, thus subdom- 
inant in (3 for high T, and ( 74 ) is the complete solution. 
Taking two derivatives one obtains /S.y{w) = —d1jR{w) 
in agreement with ( 70 ) . 



These results can be compared to simulations in the 
high-temperature region. As can be seen on figure [Tj 
they show excellent agreement. 

We thus observe that the behaviour of the model when 
random phases are added is very different in the high- 
temperature phase from what was observed in the frozen 



phase (section IV). For small random phase disorder, i.e. 



small ag, the disorder correlators stay regular at zero 
and do not develop any cusp or logarithmic singularity. 
An intuitive, physical explanation for this can be given: 
In the high-temperature phase, the fluctuations of Z are 
subdominant compared to the average Z in the m — > 
limit. Hence, even if there is a finite probability for 
fluctuations with opposite phases, a macroscopic number 
of them would need to occur simulatenously in order to 
cancel the average Z and lead to a zero of Z(w). This 
becomes infinitely improbable in the m — >■ limit. On the 
other hand, in the frozen phase, Z can be approximated 
by a two- well picture, even in the m ^ limit. Then 
there is a finite probability for the two minima to have 
opposite phases, and thus a finite probability for a shock 
where Z passes through zero. 
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VI. SUMMARY AND CONCLUSION 

In this paper, we have discussed interference effects 
in toy models of disordered systems. We considered 
one-dimcnsional models, where interference is included 
through a random complex phase on each lattice site. 

We have obtained the scaling behaviour and asymp- 
totic analytic expressions for the effective disorder cor- 
relators in the three phases of the model. For high 
temperatures small random phases do not change the 
physics, but strong random-phase disorder leads to a new 
strong-interference phase. This phase is characterized by 
a Gaussian distribution of Z centered around zero, and 
hence a finite density of zeroes of Z . For low tempera- 
tures, the system is frozen. Introducing random phases 
changes the structure of the shocks which are observed 
when a particle is "dragged" through the random poten- 
tial. There, also, zeroes of Z or, equivalently, poles of the 
complex Burgers velocity field can occur. This physical 
characterization is seen in the effective disorder correla- 
tors as a logarithmic singularity around zero. 

A few directions in which the present discussion could 
be continued come to mind. The most important aspect 
would be, certainly, to relate the physics observed here 
to higher-dimensional, more realistic models of interfer- 
ing quantum systems. In principle, one should be able to 
obtain the effective disorder correlators from field theory, 
in the frozen phase e.g. from functional renormalization 
group methods [5|. The main technical difficulty, as ap- 
parent from our toy model and a preliminary study |53| 
is the behaviour at zero: Instead of a rounding of the lin- 
ear cusp at finite temperature, we may see a logarithmic 
singularity. This makes the derivation of a field theory 
for the frozen phase in the presence of random phases a 
challenging problem. 

Another direction, which would be interesting to un- 
derstand better, is the relationship of the present results 
on the abundance of poles of the Burgers velocity pro- 
file to the pole expansion method for the solution of the 
Burgers equation |31H34| , and the pole condensation phe- 
nomena observed in 1331. 
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Note added: 

After completion of this paper we became aware of 
a very recent preprint |54j by Gredat, Dornier and Luck 
which also treats imaginary Brownian disorder motivated 
by a connection to reaction-diffusion processes. While 



the focus is different, i.e. they study the so-called Kesten 
variable which amounts to a linear potential regulariza- 
tion, while we study a quadratic well, there is agreement 
whenever the results can be compared. 



Appendix A: Moments of the partition sum with 
finite system size L and long-range correlated 
disorder 

Let us consider Q with V{x) = and 9{x) as defined 
in ( 33 ) . In this appendix, we calculate explicitely mo- 



ments of Zl. We also discuss how they can be organized 
to extract the dominant contributions for large L. 

First, we need to make some technical remarks. Con- 
sider the integral 

/,f(Ai,...,A„) := / d:ri... / " 'd:r„e^-i^--- (Al) 

k=a 1 = 1 l^j=l ^0 l=k+l l^j=k+l 



By introducing the partial sums /i^ := X^i^i -^j' the for- 



mula ( Al I can be rewritten as 



fc=0 1=0 



(A2) 



Taking a Laplace transform with respect to L, this is 
further simplified to 

/•CO 

LT{/„(A,)}(,s) = / e-^Hf(A,)dL 

n 

=n7^ 



1=0 



Mi 



(A3) 



Let us now return to the moments of ZZ*. We would 
like to evaluate 

pL pL pL pL 

{ZZ*Y = I dxi--- dxn / dyi--- d?/„ 
Jo Jo Jo Jo 

g-f (li;"j=i \xi-yi\+\yi-xj\'\xi-^j\~\vi'Vj\) 

If we assume an ordering of the 2n variables Xj and yj , 
the exponent is a linear combination of these variables. 



Hence, it is of the form of the integral (All. Each or- 



dering of the x's and y's can be mapped bijectively to a 
directed path from the lower left to the upper right cor- 
ner in a n X n lattice: The choice of an x corresponds to 
going up, the choice of a y to going right^. 



^ This generalizes straightforwardly to general moments like 
Z" (Z*)™, which give directed paths on an n X m lattice. 
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Ordering -\j -fij 
xyxy (-1,1,-1,1) (0,-l"0,-l,0) 
xyyx (-1,1,-1,1) (0,-1,0,-1,0) 
xxyy (-1,-3,3,1) (0,-1,-4,-1,0) 

xyxyxy (-1,1,-1,1,-1,1) (0,-1,0,-1,0,-1,0) 

xyxxyy (-1,1,-1,-3,3,1) (0,-1,0,-1,-4,-1,0) 

xxxyyy (-1,-3,-5,5,3,1) (0,-1,-4,-9,-4,-1,0) 

TABLE I: The Xj and m for some orderings E. 



Each ordering S can be defined by a vector with 
2n entries given by 



1 if jth variable is x 
"■J 1 if jth variable is y ' 

or equivalcntly by a vector with 2n entries 



Y^y _ ) ^ if J tti variable is y 



if jth variable is x 



(A4) 



(A5) 



Then, the resulting values of the Xj for the ordering E in 
the definition (Al) are 



-Af = (-1)^.^ 



\i=i 1=1 / 



(A6) 



A few examples for n = 2 and ?i = 3 are given in table 
|lj In order to apply (A2), we now need the partial sums 



3 = 1 3 = 1 

E {^3 - ^) 

ll=l 



4i:^3-'T.^i]-^ 

\i=i 1=1 / 



-I 2 



(A7) 



Again, see table IT 



for a few examples. 



Using formula ( A3 ) , the the Laplace transform of the 



moments can be written as: 

WW = (nir E i^-(^f) = (-0^ E ii — 



E 1=0 



(A8) 



In the interpretation of E as a directed path F = 
(wq, . . . ,W2n), with Wj on the square n x n lattice and 
Wo = (0,0), W2n — {n,n), the formula (A7| obtains a 
direct interpretation: ^fij^ is —cP, with d the distance 



, with d t 

to the diagonal. We thus obtain the interesting formula 
(setting a — 2 for simplicity in the following): 



{ZZ*y' = {n\f 



E n 



path (0,0)-i-(n,n) 



1 



(A9) 



with dw the distance of w to the diagonal. 

For the inverse Laplace transform, no closed formula 
is evident. However, from (A9| we can observe the fol- 
lowing: 

• The Laplace transform of (ZZ*) has poles at s = 

0, s — —1, s = —4, s = —9, etc. Hence, (ZZ*)"' as 
a function of L can be written as a sum of terms of 



rdcr 1, 



etc. 



• For large system sizes, the terms suppressed expo- 
nentially with L are irrelevant, and hence only the 
pole at s = needs to be discussed. 

• For each path, the pole at s = has the form -^^^ 
where z is the number of crossings of the diagonal. 
Its Laplace transform yields Hence, the dom- 
inant term for large L is given by the paths with 
the maximum number of diagonal crossings. 

• These are exactly the paths where the Xi and yi are 
paired, i.e. xyxyxyxy... or xyyxxyxy... etc. There 
are 2" such configurations. 

The final result is 



This argument provides a somewhat more detailed ex- 
planation of why the only configurations contributing to 
moments of the form (ZZ*)" are those where the replica 
are pairwise bound. When regularising the system by a 
harmonic well with mass m, we expect similar results, 
with - morally speaking - L replaced by However, 
we have not found a way to perform a more detailed 
computation using the regularization with a mass. 
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